Spin model of magnetostrictions in multiferroic Mn perovskites 
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We theoretically study origins of the ferroelectricity in the multiferroic phases of the rare-earth (R) 
Mn perovskites, i?Mn03, by constructing a realistic spin model including the spin-phonon coupling, 
which reproduces the entire experimental phase diagram in the plane of temperature and Mn-O-Mn 
bond angle for the first time. Surprisingly we reveal a significant contribution of the symmetric 
(S ■ S)-type magnetostriction to the ferroelectricity even in a spin-spiral-based multiferroic phase, 
which can be larger than the usually expected antisymmetric (S x fif)-type contribution. This 
explains well the nontrivial behavior of the electric polarization. We also predict the noncollinear 
deformation of the _E-type spin structure and a wide coexisting regime of the E and spiral states, 
which resolve several experimental puzzles. 

PACS numbers: 75.80.+q, 75.85.+t, 75.47.Lx, 75.10.Hk 



Frustrating spins in magnets often exhibit not only 
nontrivial orders but also intriguing switching and dy- 
namical phenomena associated with phase competitions. 
Effective reduction of the spin-exchange energy due to 
the frustration increases the relative importance of other 
tiny interactions such as the Dzyaloshinskii-Moriya (DM) 
interaction, single-ion anisotropy, and spin-phonon cou- 
pling. Their fine energy balance results in keen conflict of 
various states. This enables us to achieve sensitive phase 
controls and huge responses leading to new functionali- 
ties of materials, and also provides challenging issues for 
fundamental science. 

The rare-earth (R) Mn perovskites, i?MnC>3, offer one 
of the most typical examples. In this class of materi- 
als, the nearest-neighbor spin exchange is very small (~1 
mcV) relative to that in other perovskite compounds (e.g. 
~15 meV in LaTiOs) due to the cancellation of exchange 
contributions from ti g and e g orbital sectors [l|. Conse- 
quently the next-neighbor antiferromagnetic (AFM) cou- 
pling becomes comparable to the nearest-neighbor ferro- 
magnetic (FM) coupling. Their frustration gives rise to 
various competing phases including multiferroic phases 
where the frustration-induced nontrivial spin order gen- 
erates ferroelectric polarization P [2|, y] . 

Recent experiments revealed spectacular magnetoelec- 
tric (ME) phenomena in these multiferroic phases, i.e., 
magnetic- field- induced P flops [2J, |4| , colossal magneto- 
capacitance 0-Q , and electromagnons @, Q . To study 
and/or control these cross-correlation phenomena, thor- 
ough understanding of the magnetic structures, the phase 
competitions, and the coupling between magnetic and 
ferroelectric orders based on a reliable model is essential. 

However, there still remain many experimental obser- 
vations, which are not understood theoretically, 
(i) Multifurcation of the sinusoidal collinear phase at 



higher temperature (T) into four low-T phases depend- 
ing on the ionic i?-site radius (t\r) Q. 
(ii) Nontrivial r^-dependence of the magnitude and di- 
rection of P |9| . 

(iii) Apparently contradicting neutron-scattering results 
on the magnetic commensurability in the compounds 
with small tr 10H13|. 

(iv) Anomalous T-dependence of \P\ for YMnC>3 and 
ErMnOa @. 

In this Letter, we study theoretically the interplay of 
symmetric (S ■ S')-type magnetostriction (MS) and anti- 
symmetric (S x S)-type MS in a realistic spin model for 
RM11O3, and resolve all the puzzles listed above. We find 
a large (S-S) contribution to the ferroelectricity even in a 
spiral spin phase. This mechanism is generally expected 
in all the spin-spiral-based multiferroics. We start with a 
model in which the Mn S=2 spins are treated as classical 
vectors on a cubic lattice. A similar model has been ex- 
amined in Refs. [l|,|l4[ which gives the transition between 
two types of multiferroic spiral spin phases and explains 
several ME phenomena quantitatively. Here we further 
include the lattice degrees of freedom. This enables us to 
study the effect of the (S ■ S)-type MS as a source of the 
above puzzles, which has been missed thus far. 

The Hamiltonian is given by 
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where i x , i y , i z represent the integer coordinates of the 
ith Mn ion with respect to the cubic x, y and z axes. 
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FIG. 1: (a) Spin exchanges in RMn0 3 where FM (AFM) de- 
notes (anti)ferromagnetic exchange, (b) Mn(i)-0-Mn(j) bond 
in the orthorhombic lattice and local vector riij. The O ion 
is displaced from its cubic position (0) to the orthorhombic 
position (A ) at higher T. A further shift 5i,j along m,j can 
be induced at low T via the spin-lattice coupling, (c) A vs 
Jab for several R species calculated in Ref. PJ. Here A is 
normalized by the MnO bond length. 



The first term describes the spin-exchange interactions 
as shown in Fig. [TJa). The second and third terms stand 
for the single- ion anisotropy. For the local axes £j, rji 
and Ci attached to each Mn0 6 octahedron, we use the 
structural data of DyMn03 



151. The fourth term denotes 



the DM interaction. The DM vectors dij are expressed 
using five DM parameters, a a b, Pab, lab, etc and /3 C , as 
given in Ref. |16[ because of the crystal symmetry. The 
last term represents the lattice elastic term with K being 
the elastic constant. Here 5i.j is a shift of the O ion 
between ith and jth Mn ions normalized by the MnO 
bond length. Note that the O ion in the orthorhombic 
lattice is already displaced from its cubic position. We 
consider Sij as a further shift from the position at higher 
T in the presence of magnetic order at low T. 

Since the nearest-neighbor FM coupling in i?MnC>3 is 
sensitive to the Mn-O-Mn bond angle, we consider the 
Pcierls-type spin-phonon coupling J^ = J a b + J'„ h 6i. j fo r 
the in-plane Mn-O-Mn bonds where J' ab =dJ a b/dS [lj. 
We assume that the shift of O ion 8ij occurs along the 
local axis n^j directing from its cubic position (0) to the 
orthorhombic position (A ) at higher T [see Fig. Hfb)]. 
Then the positive (negative) shift decreases (increases) 
the Mn-O-Mn bond angle. 

The values of J a b, J c , Jb, D, E, and five DM param- 
eters have been microscopically determined in Ref. [1( 
for several .RM11O3 compounds. Except for Jb, they are 
nearly invariant upon the -R-site variation in the vicinity 
of the multiferroic phases. We fix J a h=— 0.8, J c =1.25, 
£=0.2, £=0.25, (a ab , p ab , 7afc)=(0.1, 0.1, 0.14), and (a c , 
/? c )=(0.42, 0.1). Here the energy unit is meV. We also 
find that very weak FM exchange J a is necessary to pro- 
duce the E phase, and adopt J =— 0.1. The value of K 
is chosen to be 500 so as to reproduce the experimental 
P in the E phase [see Fig.[3Ja)]. We obtain the value of 
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FIG. 2: Theoretical phase diagram of _RMn03 in the plane of 
Ji, and T. Here ICS denotes the incommensurate spiral phase. 
In the shaded area, the E and ICS states can coexist. Inset 
shows real-space spin configuration of the E phase. Shifts of 
the O ions due to the (S ■ S')-type magnetostriction are shown 
by gray arrows. 



J' ab from the A dependence of J a b for several R species 
[see Fig.QJc)], which gives J' ab =dJ a b/dA =2.5. 

We treat Jb as a variable which increases (decreases) 
as rn decreases (increases). This is because the ex- 
change path for Jb contains two O 2p orbitals, and the 
orthorhombic distortion, whose magnitude is controlled 
by rjj, enhances their hybridization. We find that overall 
features of the phase evolution upon the i?-site variation 
are reproduced as a function of J b even without consid- 
ering the slight independence of other parameters. 

We analyze the above model using the replica exchange 
Monte-Carlo (MC) method [18J]. Both spins and oxygen 
positions are updated in the simulation, and each ex- 
change sampling is taken after 400 standard MC steps. 
Typically, we perform 600 exchanges for a system of Af= 
48x48x6 sites along x, y and z axes with periodic bound- 
aries. We identify transition points and spin structures 
from T profiles of specific heat and spin-helicity vector 
h — jw ^2i(&i x Si+x + Si x S i+ y)/S 2 . We also calculate 
spin correlations in the momentum space by the Fourier 
transformation of spin configurations. 

We first display the theoretical T-Jb phase diagram 
in Fig. [21 which successfully reproduces the experimen- 
tally observed phase evolutions [19(. More concretely, 
the following four phases successively emerge at low T as 
Jb decreases: the A, afr-spiral, 6c-spiral, and E phases. 
In the A (E) phase, the FM (up-up-down-down) Mn- 
spin layers stack antiferromagnetically, while in the ab 
(be) spiral phase, the Mn spins rotate within the ab (be .) 
plane (Pbnm setting) to form transverse cycloids |20Ll21|. 
As T decreases, these four phases emerge with multifur- 
cation from the sinusoidal collinear phase at higher T 
where the collinear Mn spins are sinusoidally modulated 
in amplitude. The spin structure is commensurate (C) 
with q b =0.5 in the E phase, whereas it is incommensu- 
rate (IC) in the ab and be spiral phases. Importantly, 
the sinusoidal collinear state is also IC even above the E 
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FIG. 3: (Color online) (a) Polarizations vs Ji, at T— >0, i.e., (S- 
S) contribution Ps, (SxS) contribution Pas , and experimen- 
tally measured P in Eui-^Y^MnOs and Yi_ y Lu y Mn03 [jj]- 
The summation Ps+Pas reproduces the experimental P well, 
(b) Alternation of the spin angles in the ab spiral state due 
to the staggered DM vectors is illustrated in an exagger- 
ated manner where (®) denotes the positive (negative) 
c-component of the vector. Shifts of the O ions due to the 
(S ■ S)-type magnetostriction are shown by gray arrows. 



phase (e.g. (#,=0.458 for J&=2.4), and the spin-phonon 
coupling is a source of the IC-C transition with lowering 
T. 

In the ab (be) spiral phase, it has been naively be- 
lieved that the antisymmetric (S x SQ-type MS induces 
the ferroelectric polarization P\\a (P\\c) |22h24|. How- 
ever, the observed P in the ab spiral phase is much larger 
than that in the be spiral phase. For instance, the P a 
in the ab spiral phase of DyMnOa under H\\b is 2.5 
times larger than P c in the be spiral phase at P=0 [J, |5j . 
Moreover, in Euo.6Yo.4Mn03, the P a at H=0 is approx- 
imately 10 times larger than P c under iT||o, which ex- 
cludes the influence of / moments as its origin because 
of their absence [25j . This is quite puzzling since we ex- 
pect nearly identical strength of the (S x S')-type MS in 
these two phases. Recent first-principles study also sug- 
gested different mechanisms of P between the two spiral 
phases [26J. To solve this issue, we calculate the polariza- 
tion due to the (S • S)-type MS, P s = (P a , P, Pc) from 
the oxygen shifts. Because of the staggered local axes 



riij on the zigzag Mn-0 chain, P 7 (7=a, 6, c) is given by 



n 



^ = -^EK-D- 
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where (to, n) = (0, 0) for j=a, (to, n)=(l, 0) for j=b, 
and (to, n)—(i z +l, i z + l) for y=c. Here the constant n 7 
is calculated to be 4.6 xlO 5 /iC/m 2 for 7=0 and b, and 
3.3 xlO 5 /iC/m 2 for 7=c from lattice parameters using 
the point-charge model. 

In Fig. G2a), we plot calculated Ps at T— >0 as func- 
tions of Jb- Surprisingly we find a finite Ps in the ab 
spiral phase (e.g., Ps^500 fxC/m 2 for J(,=0.7), while it 



is zero in the be spiral phase. This can be understood 
as follows [see also Fig. [3jh>)] . On the in-plane chains, 
the c-axis components of the DM vectors are arranged in 
the staggered way. As a result, the spin rotation angles 
in the ab spiral become subject to an alternate modu- 
lation 27 J . Then the O ions between two spins with a 
smaller angle of </>— A0 (a larger angle of 0+A</>) shift 
negatively (positively) to strengthen (weaken) the FM 
exchange through increasing (decreasing) the Mn-O-Mn 
bond angle. These shifts generate a uniform component 
resulting in the ferroelectric polarization. In fact, the 
spin rotation angles in the be spiral are also subject to 
the alternate modulation because of the staggered a-axis 
components of the DM vectors. However the induced O 
shifts are opposite between neighboring ab planes, which 
results in their perfect cancellation. 

We also show J{,-dependence of Pas of (S x S) origin 
at T^»0 in Fig. [3ja). Since the Pas consists of two con- 
tributions, i.e. the electronic and the lattice-mediated 
ones [26J and the former one is difficult to evaluate by 
the spin model, we calculate Pas from the spin helicity 
h. Note that Pas is proportional to \h\ and the observed 
P in the be spiral phase is purely of (S x S) origin. 
In addition, we plot the experimentally measured P of 
Eui_a:Y x Mii03 and Yi-yLu^MnOa for comparison [9[, 
whose P originates purely from the Mn-spin order be- 
cause of the absence of / moments [28| . Effective r^ and 
Jb of these solid solutions are evaluated by interpolations. 

We find that the summation Ps+Pas reproduces well 
the experimental P. Here we emphasize that only the 
elastic constant K is an uncertain parameter in our 
model, and once we determine its value so as to repro- 
duce the experimental P in the E phase, the behaviors 
of P in the spiral phases are reproduced almost perfectly. 
Moreover it turns out that the (S ■ S) contribution Ps 
can be comparable to or even larger than the (S x S) 
contribution Pas in the ab spiral phase. This explains 
why P in the ab spiral phase is much larger than that in 
the be spiral phase. 

Next we discuss the E phase. Interestingly we find a 
finite c component of the spin-helicity vector h in this 
phase, indicating that its spin structure is not collinear 
in contrast to what we have believed so far, but its up- 
up-down-down structure is subject to a spiral modula- 
tion within the ab plane. The inset of Fig. [2] illustrates 
the real-space spin configuration of the E phase, which 
indeed shows the elliptically deformed a6-planc cycloid. 
Calculating the T-dcpcndcnce of the expectation value 
for each term in the Hamiltonian, we find that the single- 
ion anisotropy or alternation of the in-plane easy mag- 
netization axes due to the d 3x 2_ r 2 / d 3y 2_ r 2-ty i pe orbital 
ordering is an origin of the cycloidal deformation. This 
predicted deformation should be confirmed in a future 
polarized neutron-scattering experiment. 

With dominant up-up-down-down spin 6-axis compo- 
nents, the O ions between nearly (anti)parallcl Mn-spin 
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FIG. 4: (Color online) (a) Calculated spin-correlation func- 
tions in the momentum space when the E and incommen- 
surate spiral (ICS) states coexist. Here S 1 (k) denotes the 
correlation function for the spin 7-axis components, (b) Cal- 
culated T dependence of Ps in the pure E state and that in 
the ICS state for 2.6< feT <5.1. Since the Ps in the coexist- 
ing regime necessarily becomes smaller than that for the pure 
E state, an anomaly can be observed in the T profile of P as 
indicated by the dashed line. 



pairs shift negatively (positively) to modulate the FM 
exchanges, which results in the ferroelectric polarization 
(see inset of Fig. [5]) |29j, |30j. In Fig. [2ta), we indeed see 
a very large Ps (^4600 /iC/m 2 ) in the E phase. We also 
expect a small but finite (S x S) contribution Pas due to 
the cycloidal deformation. The calculated Ps, Pas an d 
their sum Ps+Pas in the E phase are invariant upon the 
Jb variation in agreement with the experiment [9j . 

Finally we discuss the coexistence of the E state and 
the IC-spiral (ICS) state. In the shaded area of Fig [51 
we obtain coexisting solutions when we perform the MC 
calculation starting from the ICS state as an initial con- 
figuration. The spin-phonon coupling or the (S ■ S^-type 
MS make the transition between ICS and E phases of 
strong first order. As a result, although the energy com- 
parison gives the transition line between them as indi- 
cated by the solid line, the ICS state has a deep energy 
minimum even in the E phase. This can result in their 
coexistence. Such a coexistence easily occurs in reality 
since the system enters into the E phase necessarily via 
the IC sinusoidal collincar phase with lowering T. 

To understand the contradicting neutron-scattering re- 
10Hl3|. we calculate the spin-correlation functions 



suits 

in the momentum space for a coexisting solution obtained 
in the MC calculation. We find only the spin &-axis com- 
ponent has sharp peaks at ^=±0.458, while the other 
components have very small peaks as shown in Fig. SJa) . 
This seems as if the spin structure were IC collinear, 
which leads to the confusion. Observations of the IC wave 



numbers q b ~0.43 in PMn0 3 with P=Ho [llj and Er 13 1 
as well as a report of the IC collinear state in YMn03 12 1 
can be attributed to this coexisting state, while a report 
of the commensurate ^=0.5 in H0M11O3 [lfj is to the 
pure E state. 

When 1.8< Jb <2.5, the energy minimum of the ICS 
state disappears as T decreases. In this case, an anomaly 
should appears in the T profile of P or dielectric constant 



[see Fig.Ufb)]. A recent experiment indeed observed such 
an anomaly in YMn03 and ErMnOa [9| , which strongly 
evidences the coexistence. The coexistence together with 
the (S ■ S) contribution in the E phase should be seri- 
ously considered also when we interpret the experimental 
results for PM11O3 with R=Y, Ho, ...,Lu, like the strange 
electromagnon spectra in the optical spectroscopy 31 1 . 

In summary, we have theoretically studied the mag- 
netic structures and the ME coupling in PMn03 by using 
a realistic spin model including the spin-phonon coupling. 
We have succeeded in reproducing the entire phase dia- 
gram of PMnOs for the first time, and have revealed the 
cooperative contributions of symmetric (S ■ S)-type and 
antisymmetric (S x S)-type MSs to the ferroelectricity 
in the ab spiral phase. This mechanism is generic and 
is relevant to all the spin-spiral multifcrroics. We have 
also found the cycloidal spin deformation in the E phase, 
and the coexistence of the E and ICS states. On these 
basis, the nontrivial behavior of P and several puzzles in 
the experiments have been explained. Our model gives 
a firm basis for studying and controlling the intriguing 
cross-correlation phenomena in PMn03. 
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